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All biological processes are controlled by complex systems of enzymatic chemical reactions. 
Although the majority of enzymatic networks have very elaborate structures, there are many 
experimental observations indicating that some turnover rates still follow a simple Michaelis- 
Menten relation with a hyperbolic dependence on a substrate concentration. The original 
Michaelis-Menten mechanism has been derived as a steady-state approximation for a single- 
pathway enzymatic chain. The validity of this mechanism for many complex enzymatic 
systems is surprising. To determine general conditions when this relation might be observed 
in experiments, enzymatic networks consisting of coupled parallel pathways are investigated 
theoretically. It is found that the Michaelis-Menten equation is satisfied for specific relations 
between chemical rates, and it also corresponds to the situation with no fluxes between 
parallel pathways. Our results are illustrated for simple models. The importance of the 
Michaelis-Menten relationship and derived criteria for single-molecule experimental studies 
of enzymatic processes are discussed. 
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I. INTRODUCTION 

It is known that all chemical reactions in biological systems are catalyzed by protein enzymatic 
molecules.— 1 " 2 Fundamental understanding of cellular processes cannot be accomplished without 
determining how involved enzymatic networks function. First catalytic mechanism that involves 
enzyme molecules has been proposed almost a century ago by Michaelis and Menten,— and since 
it became one of the most used and celebrated relations in biochemical and biophysical studies 
of natural phenomena. Similar relations between chemical rates and substrate concentrations 
have been observed in many enzymatic systems.— ~— However, why the Michaelis-Menten (MM) 
mechanism, i.e., the hyperbolic dependence on the concentration of substrate, is working even for 
some very complex enzymatic networks is still not well understood. 

The original MM mechanism has been derived as a steady-state approximation to a simple 
single-pathway enzymatic process with irreversible creation of product molecules as shown in a 
scheme at Fig. 14 In addition, it was assumed that the initial concentration of enzyme was much 
smaller than the initial concentration of the substrate, and that not too many product molecules 
have been produced.— Then the rate of the catalyzed reaction can be written in terms of the rate 
constants shown in Fig. 1 as 

c + K M w 

where c is the concentration of substrate molecules, Km = (^2 + k_i)/k\ is known as the Michaelis 
constant and the rate ki is called the catalytic rate. The rate of enzymatic reaction in the MM 
mechanism is proportional to the concentration of the substrate molecules at low concentrations 
while it saturates at high concentrations. Surprisingly, similar behavior is also observed experi- 
mentally in many enzymatic systems with much more complex topology of chemical transitions 
than one shown in Fig. 1.— — Several theoretical studies have been presented in order to understand 
why turnover rates in complex enzymatic networks might follow the simple MM relation.— _ — It 
has been argued that the MM equation will be satisfied for fluctuating enzymatic systems when 
conformational transition rates are very large or very low in comparison with catalytic rates, al- 
though with redefined definitions of catalytic rates and Km constants.— 1 ^ A general approach to 
analyze enzymatic kinetics, based on the flux balance method, has been developed recently— and 
it suggests that the MM relation will not be observed when the detailed balance is broken. How- 
ever, in all current approaches only enzymatic networks with several irreversible transitions have 
been considered. In addition, general expressions for the rate of product formation have not been 
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derived, and explicit conditions on validity of the MM equation have not been obtained. 

In this paper, we present an alternative theoretical method of analyzing complex enzymatic 
networks based on solving discrete-state stochastic master equations. This approach has been 
utilized successfully to describe dynamic properties and mechanisms of motor proteins, and also in 
the analysis of other single-molecule experiments.— - — Specifically, we consider enzymatic networks 
consisting of coupled parallel pathways with all chemical transitions assumed to be reversible. Our 
method allows us to describe explicitly conditions when the MM relation for the turnover rate 
is satisfied. These conditions are analogous to the requirement that the flux between parallel 
pathways vanishes at each state. Using simple models to explain our method we argue that all 
turnover rates behave qualitatively similar to the MM relation. 

The paper is organized as follows. In Sec. II our theoretical method is presented for general 
coupled parallel enzymatic networks. Examples to illustrate obtained results are given in Sec. III. 
Finally, Sec. IV provides a summary of our findings. 

II. THEORETICAL ANALYSIS 

As an example of complex enzymatic networks a system made of several coupled parallel bio- 
chemical pathways, as shown in Fig. 2, will be analyzed. We adopt here a single-molecule approach, 
and to simplify calculations all derivations will be made for the network of two coupled parallel 
pathways (see Fig. 2), although the analysis can be easily extended to many coupled parallel path- 
ways. This situation corresponds to a single enzymatic molecule that can be found in one of two 
conformations, and it catalyzes the reaction of the substrate transformation to the corresponding 
product in both conformations but with different rates. 

It is assumed that in each pathway there are N discrete chemical states per each enzymatic 
cycle. The enzyme molecule can be found in the state i from which it can move forward or backward 
along the same reaction channel with rates u% or Wi, respectively, or it can change its conformation 
by moving to the state i in the second pathway with the rate 7$: see Fig. 2. Similarly, if the 
molecule is in the state i of the second reaction channel it can make forward (backward) transition 
with the rate oti (/%), or it can change its conformational state by going to the first pathway with 
the rate <5j, as illustrated in Fig. 2. We also assume that the enzyme molecule binds to the substrate 
in the states i = 0, i.e., the corresponding transition rates are proportional to the concentration 
of substrate molecules c, uq = k\c and ao = k^c. Let us define a function P^ k \t) as a probability 
to find the enzyme molecule in the state i at the reaction channel k (k = 1 or 2) at time t. The 



temporal evolution of this enzymatic network can be described by a system of master equations, 

= u^Mt) + w i+1 p£\(t) + Siifit) - ( Ui + w i + 7i )^ (1) (t); (2) 

= a^P^t) + P W P%\(t) + 7 ^ {1) (t) - (a t + ft + ft)jf >(t). (3) 

Because of periodicity, it is convenient to define a new function, 

RW(t) = J2pM N (t), (4) 
i 

with / being any integer number. Then the conformational transition flux between the states i 
(i = 0, 1, . . . , N — 1) in two pathways can be written as 

J i = li Rf ) -5iRf\ (5) 

When the system reaches a stationary state (t — > oo), Master equations ([2]) and ([3]) transform 
into the following expressions, 

= Ul ^Rf\ + Wi+rR^t) -Ji- fa + Wi)n!p; (6) 

= ai^Rf\ + fa+iR^ + Ji- {on + Pi)Rf ] . (7) 
In addition, at large times total conformational flux between states in two pathways must disappear, 

jV-l 



£ J * = °- ( 8 ) 



i=0 

Now we can define new effective rates, 



and 



Ui = Ui + aiY, (9) 



Wi=Wi+f3i — . (10) 



Using these definitions and the expression for the flux ([5]) one can obtain from Eq. ([6]) 

= CTii^i - mR^ - Km - UiR^ + m^iR^l + Ki, (11) 



where a new function K; is defined as 



Ki = -iJi - ^Ji+i- (12) 

Oi Oi+l 



Since Eq. (fTT]) must be satisfied for any chemical state i it leads to 



where C is some unknown constant. Using periodicity Eq. (113j) can be easily solved, yield 

(i) _ Cn + ri(K) 



R 



1 - r 



with 



AT-1 — 



«=0 



, , AT-1 k 



» k=l 1 = 1 J 



and 



JV-l fc 



fc=i i=i 

The unknown constant C can be found from the normalization condition, namely, 

N-l 



E{ R i+ R i) = 1 



which produces 

^_ (i + L)(i-r)-5(K) 
c - S ' 

where new auxiliary functions are defined as 

1=0 



i=0 
and 

2V-1 



N-l 

6; 



i=0 

Substituting Eq. ([19]) into Eq. (fb^) it can be shown that 

(1) _ /i + L _ _m\ r . n(K) 



R R(1-T)J (i-r) 
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These expression for the probability to be found in the state % of the first reaction channel at 
steady-state conditions allow us to compute all dynamic properties of the system. The turnover 
rate for this enzymatic network, 

N-l N-l 

v = j2(-"i- + E fa - A)*i 2) . ( 24 ) 

which can be written as 



i=0 i=0 



N-l N-l 

?(!) 



4=0 1=0 



Then it can be shown that 



JV-l 



^2(ui-m)ri = N(l-T), (26) 



and 



i=o 



N-l N-l 



Y J {ui-w i )r l {K) = {l-T)Y J K l . (27) 



i=0 i=0 

It leads to a compact expression for the chemical reaction rate, 



v= N^ + Na-r)L + NRp m 

o o o 



We can repeat derivations if we choose another set of effective rates, 



and 



5' 

&i = at + Ui — , (29) 

H 



(3 i = f3 l + wA. (30) 



Comparing them with Eqs. Q and (fit))) we conclude that 



ok ^ $i 



(31) 



Again the equation for the turnover rate can be derived with this set of effective rates. Combining 
both expression, we arrive to the final formula for the reaction speed of the enzymatic process 
presented in Fig. 2, 
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with 

N-l k 



(M) = l\M i+ J2Mi + kf[^), (33) 



Ui 



N-l , , 

w = j i+? r<(n (34) 

and functions Mj are defined as 

\<>i 7i/ V^+l H+lJ 

Eq. (|32p is an exact expression for the turnover rate for general enzymatic network made 
of coupled parallel pathways. It depends on chemical rates for individual transitions as well as 
conformational fluxes Jj. It should be noted that these fluxes are not independent parameters. 
Explicit forms of conformational fluxes can be found utilizing Eqs. © and also from the balance 
of total fluxes through each chemical state,— 

Ui-iflS - WiR? = Ji + mR? - w i+1 R^ v (36) 

c^iifi - ftijp } + J l = ai Rf ] - fa+iR^ • (37) 

It is convenient to analyze the expression (|32[) for the overall chemical rate of the coupled 
enzymatic network if it follows or not the MM dependence. It can be easily shown that the first 
term has a hyperbolic dependence on the concentration of the substrate and it does not depend on 
conformational fluxes. At the same time, the second and the third terms do depend on Jj and they 
are not satisfying the MM relation for any non-zero conformational fluxes. Thus the condition for 
the turnover rate to follow the MM dependence is when last two terms in Eq. (I32p vanish, which 
can only take place when conformational fluxes disappear. This condition can be written explicitly 
in terms of relations between chemical rates, 

a-jr^ (38) 

or in the different form as 

Uj, _ Cti Wi _ Pi 
7t °i 7i Oi 

This is the main result of our work. A general coupled parallel enzymatic network will follow the 
MM relation when rates of coupled pathways satisfy the equation (I38h for each chemical state. This 
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is a clear and explicit criterion on the MM dependence of the enzymatic cycle in such complex 
systems. 

All previously discussed situations of the MM behavior for the turnover rate in enzymatic 
networks ^ 101 " 12 can be described by our general approach. It has been argued that conforma- 
tionally fluctuating enzymes still obey the MM relationship for 1) quasi-static conditions (when 
conformational rates are very slow in comparison with other chemical transitions), and for 2) 
quasi-equilibrium conditions (when conformational rates are very fast). The last case has been 
also analyzed in Ref. 16 . The first situation corresponds to 7$ <C 1 and Si <C 1. Then Eqs. ([39]) 
are valid for any value of other chemical rates (as long as they are larger than conformational 
transitions rates). The second case describes very large 7^ and Si, and again Eqs. (i39j) predict that 
it might happen for any value of chemical rates Ui, Wi, a.i and /3% . It was also found 1 - that the MM 
relation still holds when only one group of the conformationally-related states i is found in quasi- 
static or quasi-equilibrium conditions. Since N = 2 periodic system has been used for analyzing 
fluctuating enzymes, then from Eq. (JHJ) one can easily conclude that there are two conformational 
fluxes, J\ and J%, and they are related via J\ + J2 = 0. If one of them goes to zero the second flux 
must also vanish, and the MM relationship is recovered. Cao has argued that the MM dependence 
is observed when the detailed balance is not broken.— For the system of coupled parallel pathways 
it means that the overall circular current in each loop (see Fig. 2) is zero. The ratio of clockwise 
and counterclockwise currents for any loop can be written with the help of Eqs. (|38p as 

ugj+iPi+iSj _ Tj+iA+i _ 1 ^ 
aiS i+1 w i+1 ^i S i+ iw i+ i 

which indicates that the loop current vanishes when the MM relation is observed. 

It is important to note also that criteria given in Eqs. (]38p suggest that the MM dependence 
is found for systems where each state has similar free-energy landscape near each chemical state. 
Then the enzymatic network of coupled parallel reaction channels can be effectively viewed as a 
single enzymatic pathway with properly rescaled transitions rates which by definition follow the 
MM behavior. 

III. EXAMPLES 

To illustrate our theoretical method let us present some explicit results for the N = 2 system 
which is the most relevant for analyzing fluctuating enzyme systems.— The dynamic proper- 
ties of enzyme molecule in such system has been already obtained.— The general formula for the 
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turnover rate is given by**- 

V = (l/n){kicux - w wi)[(k 2 c+^ )6 1 + (a 1 +Pi)5o + 6 6 1 ] (41) 
+ (k 2 cai - /3 /3i)pic + wohi + (tii + wi)lo + 7o7i] 
+{k x ca\ - wo/3i)5 7i + (k 2 cu x - Aru7)7oM> 

where the function is defined as 

n = (feic + to + 'yD)[(fc2c + A))*i + (Q!i + /9i)7i] ( 42 ) 
+(ui + wi + 7i)[(fe 2 c+ /3o)7o + ("i + /3i)<5 ] 
+(A: 2 c + /3 + 6 )[(kic + tu )7l + C"i + 
+(«! + /3i + (5i)[(/cic + wo)£o + (ui + ^1)70]. 

The conformational flux is equal to 

J = -J 1 = J= (l/fi)[(ui + W!)(k 2 c + Pohoh ~ (kic + wdiax+PdrfxSo]. (43) 

When the conformational fluxes are zero the turnover rate has the MM dependence, 

y (Mo + fc27o)( n i^i + «i7i)c - (wq^o + /3q7q)(wi<5i + ^171) 

[(kic + w )S + (k 2 c + ^0)70] (71 + <5i) + [«i +^1)^1 + (ax + #1)71] (70 + ^0)' 

Reaction rates for different sets of parameters are presented in Fig. 3. The dependence of all 
enzymatic rates on the substrate concentration is qualitatively similar: increasing function at small 
c and saturation to a constant value for c> 1. There is nothing special in the MM relation (black 
curve in Fig. 3). This can be easily understood because the turnover rate is equal to the ratio 
of polynomial functions of c. For the MM case both numerator and denominator are just linear 
functions. At large concentrations of the substrate molecules transitions from the states i = 1 
become rate-limiting leading to effectively concentration-independent behavior. One could also 
observe that changing transition rates for i = does not affect much overall enzymatic rates in 
comparison with the MM relation (green and blue curves in Fig. 3). However, modifying transition 
rates associated with i = 1 can significantly alter reaction rates (red and brown curves in Fig. 3). 
In addition, analyzing all reaction rates curves we conclude that many experimental observations 
for systems that could analyzed by iV = 2 coupled periodic models might be easily assumed to 
follow effectively the MM relation, especially if experimental errors are taken into account. 

The conformational fluxes for different sets of parameters are shown in Fig. 4. In the MM 
case the J = as expected, while in other cases the complex behavior is observed. At large 
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concentrations all conformational fluxes tend to zero because in this case transitions from i = 1 
are becoming rate limiting, and they do depend on c. In some situations the magnitude of the 
conformational flux reaches a maximum for the specific concentration of the substrate molecules. 

IV. SUMMARY AND CONCLUSIONS 

A new theoretical method of analyzing complex enzymatic networks is developed. It is based 
on solving explicitly master equations for discrete-state stochastic models. Applying this approach 
for systems made of coupled parallel enzymatic pathways, we derived explicit criteria on when 
the MM behavior might be observed. It leads to special relations between all chemical rates. 
This method allows to explain all previous theoretical observations on the validity of the MM 
relationship in complex enzymatic systems. Theoretical analysis is illustrated for simple models 
relevant for conformationally fluctuating enzyme molecules. It is found that general behavior of 
enzymatic rates is qualitatively similar to the MM equation. Our theoretical method argues that 
turnover rates can be viewed as a ratio of polynomial functions of the substrate concentrations. 
The MM relationship is observed when these polynomials are linear functions of concentrations. 

The presented theoretical method might be a powerful tool for analyzing single-molecule exper- 
iments since it allows to compute all dynamic properties for complex enzymatic networks. It will 
be interesting to test experimentally developed criteria for coupled enzymatic pathways. It will be 
also important to extend this approach to other enzymatic systems and for computations of other 
experimentally observed quantities such as dwell-time distributions and diffusion constants. 
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Figure Captions: 

Fig. 1. A reaction scheme for the simple enzymatic process that can be described by the Michael- 
Menten (MM) mechanism. An enzyme molecule E reacts reversibly with a substrate molecule S 
to produce an intermediate molecule X which can irreversibly transition to a product molecule P 
and the enzyme E with corresponding rates. 

Fig. 2. Enzymatic network with two coupled parallel pathways. Both pathways have N discrete 
states per each enzymatic cycle. The enzyme molecule in the state i can transition forward (back- 
ward) with the rate Uj (u>j) if found in the pathway 1, while in the second channel the forward 
(backward) rates are a, (/%) (with i = 0, 1, . . . , N — 1). The conformational transitions between 
two i-th. states in different pathways are given by rates 7« and <5j. 

Fig. 3. The turnover reaction rate as a function of the concentration of substrate molecules for 
N = 2 periodic coupled parallel enzymatic network. Calculations have been performed utilizing 
Eq. (|41|) for the following parameters for all curves: k\ = 10 /xM -1 s , &2 = 1 /uM _1 s -1 , u\ = 5 
s -1 , wi = 70 = 1 s -1 , /3q = Pi = 0.1 s _1 and 5\ = 0.5 s . In addition, for the brown curve we 
used a.\ = 5 s™ 1 , 5q = 0.1 s _1 , 71 = 5 s _1 and wq = 1 s _1 ; for the green curve we used a\ = 0.5 
s , 5q = 1 s -1 , 71 = 5 s _1 and wq = 1 s _1 ; for the black curve we used a\ = 0.5 s _1 , 5q = 0.1 s _1 , 
71 = 5 s" 1 and wo = 1 s _1 ; for the blue curve we used a\ = 0.5 s _1 , 6q = 0.1 s _1 , 71 = 5 s -1 and 
wq = 10 s ; and for the red curve we used a\ = 0.5 s -1 , 5q = 0.1 s _1 , 71 = 50 s -1 and Wq = 1 
s . The Michaelis-Menten case is described by the black curve. 

Fig. 4. The conformational flux as a function of the concentration of substrate molecules for N = 2 
periodic coupled parallel enzymatic network. Calculations have been performed utilizing Eq. (|43p 
for the following parameters for all curves: k\ = 10 /iM _1 s _1 , k2 = 1 [iMT 1 s _1 , u\ = 5 s _1 , 
w i = 7o = 1 s_1 j Po = Pi = 0.1 s^ 1 and <5i = 0.5 s" 1 . In addition, for the brown curve we used 
Oi\ = 5 s _1 , 5o = 0.1 s _1 , 71 = 5 s _1 and Wq = 1 s _1 ; for the green curve we used «i = 0.5 s _1 , 
5o = 1 s _1 , 71 = 5 s _1 and wq = 1 s _1 ; for the blue curve we used a\ = 0.5 s -1 , <5o = 0.1 s _1 , 
71 = 5 s -1 and wq = 10 s^ 1 ; and for the red curve we used a± = 0.5 s _1 , 5q = 0.1 s _1 , 71 = 50 s _1 
and wo = 1 s _1 . For parameters describing the MM relationship the conformational flux is equal 
to zero. 



12 




Figure 1. Kolomeisky 




Figure 2. Kolomeisky 




Figure 3. Kolomeisky 




Figure 4. Kolomeisky 



